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Abstract 

We calculated the vector, axial-vector, scalar and pseudo-scalar two-point 
functions up to two-loop level in the low-energy effective field theory for three 
different QCD-like theories. In addition we also calculated the pseudo-scalar 
decay constant Gm- The QCD-like theories we used are those with fermions in 
a complex, real or pseudo-real representation with in general n flavours. These 
case correspond to global symmetry breaking pattern of SU{n)L x SU{n)ii — 
SU{n)v, SU{2n) SO{2n) or SU{2n) Sp{2n). We also estimated the S 
parameter for those different theories. 
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1 Introduction 

The different global symmetry breaking patterns of QCD-like theories with a vector- 
like gauge group have been summarized in [H |2l |3] around 30 years ago. The global 
symmetry and its spontaneous breaking depend on whether the fermions live in a 
complex, real and pseudo-real representation of the gauge group. For n identical 
fermions this corresponds to the symmetry breaking pattern SU{n)L x SU{n)fi — t- 
SU{n)y, SU{2n) — )■ S0{2n) and SU{2n) — > Sp{2n) respectively. These theories can 
be used to characterize some of technicolor models with vector-like gauge bosons. 
QCD-like theories are also important in the theory of finite baryon density. Here the 
real and pseudo-real case allow to investigate the mechanism of diquark condensate 
and finite density without the sign problem. A main nonperturbative tool in study- 
ing strongly interacting theories is lattice gauge theory. Numerical calculations are 
performed at finite fermion mass and need in general to be extrapolated to the zero 
mass limit. In the case of QCD Chiral Perturbation Theory (ChPT) is used to help 
with this extrapolation. Our work has the intention of providing similar formulas 
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for the QCD-like theories using the effective field theory (EFT) appropriate for the 
alternative global symmetry patterns. 

These EFT have been used at lowest order (LO) [1] with earlier work to be found 
in [5l El [7] and some studies at next-to-leading order (NLO) have also appeared 
[HI m ID] • The former two are the usual QCD case with n flavours. In our earlier 
papers [HI [T2] we have systematically studied the effective field theory of these three 
different QCD-like theories to next-to- next-to-leading order (NNLO). We managed to 
write the EFT of these cases in an extremely similar form. We calculated the quark- 
antiquark condensates, the mass and decay constant of the pseudo-Goldstone bosons 
and meson-meson scattering [12]. In this paper we extend the analysis to two- 
point correlation functions. We obtain expressions for the vector, axial-vector, scalar 
and pseudo-scalar two-point functions as well as the pion pseudo-scalar coupling Gm 
to NNLOQ or order 

In our earlier work [HI [12], we called the three different cases QCD or complex, 
adjoint or real and two-colour or pseudo-real. In this paper we use only the latter, 
more general, terminology. 

One motivation for this set of work was the study of strongly interacting Higgs 
sectors, reviews are [131 [H]. For any model beyond the Standard Model, passing the 
test of oblique corrections, or precision LEP observables, is crucial [T5l [T6] . Over the 
years, the impact of the oblique corrections in those models have been studied quite 
intensively but in strongly interacting cases mainly an analogy with QCD has been 
invoked. Lattice gauge theory methods allow to study strongly interacting models 
from first principles. The contributions from these theories to the S'-parameter can be 
calculated using the two-point functions studied here and our formulas are useful to 
perform the extrapolation to the massless case. This was in fact the major motivation 
for the present work but we included the other two-point functions for completeness. 

The paper is organized as follows. In Section |2] we give a brief introduction to 
EFT for the three different cases. Section |3] is the main part of the paper. We define 
the fermion currents and the two-point functions in Section 13.11 In Sections 13.21 
to 13. 5[ we present the calculation of vector, axial-vector, scalar, pseudo-scalar two- 
point functions. In Section |H we discuss the oblique corrections and the S-parameter. 
Section El summarizes our main results and we present the definition. 

2 Effective Field Theory 

In this section we briefly review the EFT of QCD-like theories, the details can be 
found in the earlier paper [11]. The basic methods are those of Chiral Perturbation 
Theory [T71 HE]. The counting of orders is in all cases the same as in ChPT, we count 
momenta as order p and the fermion mass m as order . 

2.1 Complex representation: QCD and CHPT 

The case of n fermions in a complex representation is essentially like QCD. The 
Lagrangian with external left and right vector, scalar and pseudos-calar external 

^We use LO, NLO and NNLO as synomyms for order p^, order and order calculations even 
if the order vanishes. 
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sources, l^,r^,s and p, is 

^ = Qu^YDf^qii + qmiYDf^qm + quI^^rnQLj + qmYrt^ijqRj 

-QmMijqij -QuMljqRj i,j = l,2,...,n. (1) 

The covariant derivative is given by D^g = — iG^q, and the mass matrix M. = 
s — ip. The sums shown are over the flavour index. The sums over gauge group 
indices are imphcit. 

The Lagrangian ([T]) has a symmetry SU{n)L x SU{n)fi which is made local by 
the external sources [SI [IB]- The quark- ant i- quark condensate (qq) breaks SU{n)L x 
SU{n)R spontaneously to the diagonal subgroup SU{n)v. According to the Nambu- 
Goldstone theorem, — 1 Goldstone Bosons will thus be generated. We add a small 
fermion mass m explicitly by setting s = m + s. This mass term explicitly breaks 
the symmetry SU{n)L x SU{n)ji down to SU{n)y as well and gives the Goldstone 
bosons a small mass. 

The Goldstone boson manifold SUiji)^ x SUin)^/ SU{n)v can be parametrized 

by 

M = exp (^-^tt'^T'^^ a = l,2,...,n2-l. (2) 

The T"^ are the generators of SU{n) normalized to {T"'T^) = 6"'^. The notation (A) 
stands for the trace over flavour indices, u transforms under Ql ^ Qr ^ SU{n)L x 
SU {71)11 as u — )■ gRuK^ = hug\ where h is the "compensator" and is a function of 
ql and qr. The methods are those of [1^, but we use the notation of [201 |2T]. We can 
construct quantities which transform under the unbroken group if as : O — )■ hOK^ 



u 



V^O = d^O + V^O-OV^, 

X± = u^XU^ ± ux^u , 
f±^^ = ul^,yv) ± v)r^yU . (3) 

The field strengths 1^^ and r^^ are 

/ = (j 1 - (j 1 ~ 711 /I 
I'fj.u '^^I'U '^u''iM '[''/^) I'Ui ; 

rf,u = d^r^ - d^r^ -i[r^,, r^] . (4) 
The covariant derivative contains 

= ^[u^d^ - ir^)u + u{d^ - lf,)u^] ■ (5) 

X contains the matrix Ai, which is the combination of scalar and pseudo-scalar 
sources 

X = 2BoM=2Bois-ip). (6) 

Using the quantities in ([3]), we can find the leading order, j9^, Lagrangian which is 
invariant under Lorentz and chiral symmetry: 

/:2 = ^(mX + x+)- (7) 

The subscript "2" stands for the order of p"^. The p'^ and p^ Lagrangian will be 



explained in Section 1231 
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2.2 Real and Pseudo-Real representation 

The case of n fermions in a real or pseudo-real representation of the gauge group 
we can treat in a similar way as the complex case. In the real case, the global 
symmetry breaking pattern is SU{2n) — )• S0{2n), and the number of generated 
Goldstone bosons is 2n^ + n — 1. In the pseudo-real case, the symmetry breaking is 
SU{2n) — > Sp{2n), and the number of generated Goldstone is 2n^ — n — 1. In both 
cases anti-fermions are in the same representation of the gauge group and can be put 
together in a 2n vector q, see |lT] for more details. 

The condensate can now be a diquark condensate as well as a quark- ant iquark 
condensate. Our choice of vacuum corresponds to a quark- ant i- quark condensate. In 
terms of the 2n fermion vector q they can be written as 

Real : {q^CJsq) + h.c. Js = 

Pseudo - Real : {qa^apC JAq/s) + h.c. ^ I o" J ' 

Here C is the charge conjugation operator. Js and Ja are symmetric and anti- 
symmetric 2n X 2n matrices, / is the nxn unit matrix. Since Js and Ja often appear 
in the same way in the expressions, we use J for both cases unless a distinction is 
needed. 

The generators, T", of the global symmetry group SU{2n) can be separated 
into belonging to the broken, or unbroken part, Q". They satisfy the following 
relations with J: 

JQ" = -Q^^J, JX" = X^^J , (10) 
The Goldstone boson manifold can be parametrized with 

U = uJu^ gUg'^ , with u = exp • (11) 

where J = Js and a runs from 1 to 2n^ + — 1 for the real case and J = Ja and a 
runs from 1 to 2n^ — n — 1 for the pseudo- real case. 

In our earlier paper [11], we constructed quantities similar to those in (jSHS]) 

= i[u\d^ - iVf,)u - u{d^ + iJV^ J)u^] , 

= ^[u\d^-tV^)u + u{d^ + tJV^J)u^]. 

f±iMu = JuV^^uU ±uVf,^u^ , 

u^xJu'^ i uJx^u . (12) 



X± 



The 2n x 2n matrix includes the left and right-handed external sources 

V. - ( -,r ) (13) 

and V^y is the field strength 

V = d^y, - d,V^ - % {Vy, - V,V^) . (14) 
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X include the matrix Ai via x = 2BqM. [H]. Those quantities behave similarly as 
those ([3]) if we take 

-JV^J^l^^ V,^r,. (15) 

With this correspondence, the Lagrangian of the real and pseudo-real case has the 
same form as the complex one. However one has to remember there are differences 
in the generators, external sources, coupling constants,. . . . Anyway, now we can use 
the techniques of ChPT to perform the calculations. 

2.3 High Order Lagrangians and Renormalization 

Using Lorentz and chiral invariance, we can write down the EFT lagrangian [18] 
for all three cases using the quantities listed in ([3]) and f|T2l) : 



+L4(m%)(x+) + L5(m%x+) + L^{X+? + Lj{x-? + \ls{xI + X-) 
-iU{U^,u>^u^) + \L,,{fl - f_) + H,{l^,F^ + r^.r'^O + H^ixX^) . (16) 

To do the renormalization, we use the ChPT MS scheme with dimensional regular- 
ization [181 El EI] . The bare coupling constants Li are defined as 

L, = (c/x)'-^[r,A + L[(/i)] , (17) 

where the dimension d = 4 — 2e, and 

^ = 167r2(d-4) ' ^^^^ 
Inc = -^[ln47r + r'(l) + l] . (19) 

The coefficients Fj for the complex case have been obtained in [8j, for the real and 
pseudo-real case we have calculated them earlier in [TT] . However, there are mistakes 
in the coefficients of Lq,Liq and Hi in the Table 1 of [TT]. These mistakes had no 
effects on our previous calculations. We therefore list all the coefficients here again 
in Table [T] 

The Lagrangian for the complex case and general n has been obtained in [20] , 
it contains 112+3 terms. The divergence structure of the bare coupling constants Ki 
in the p^ can be written a^ 



2(d-4) 



IGvr^ 



(20) 



The coefficients Tf'\ F^^^ and F^-^'' for the complex case have been obtained in [21] . 

For the real and pseudo-real case, the p^ Lagrangian has the same form as in the 
complex case but some terms might be redundant. The divergence structure as given 
in (120 p still holds but the coefficients are not known. One check on our results that 
remains is that all the non-local divergences cancel. 



^The Ki have been made dimensionless by including a factor of explicitly in the order 
Lagrangian. 
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i 


complex 


real 


pseudo-real 





ra/48 


(n + 4)/48 


[n - 4) /48 


1 


1/16 


1/32 


1/32 


2 


1/8 


1/16 


1/16 


3 


n/24 


(n-2)/24 


(n + 2)/24 


4 


1/8 


1/16 


1/16 


5 


n/8 


n/8 


n/8 


6 


(n^ + 2)/(16n2) 


(n2 + l)/(32n2) 


(n^ + l)/(32n2) 


7 











8 


(n^ -4)/(16n) 


+ n - 2)/(16n) 


(n2 - n - 2)/(16n) 


9 


n/l2 


(n + l)/12 


(n - 1)/12 


10 


-n/12 


-(n+ 1)/12 


-(n - 1)/12 


1' 


-n/24 


-(n+ l)/24 


-(n - l)/24 


2' 


(n2 -4)/(8n) 


+ n - 2)/(8n) 


(n^-n -2)/(8n) 



Table 1: The coefficients Fj for the three cases that are needed to absorb the diver- 
gences at NLO. The last two lines correspond to the terms with Hi and H2. This is 
the same as Table 1 in [11] but with the error for Lg, Liq and Hi corrected. 

3 Two-Point Functions 
3.1 Definition 

The effective action of the fermion level theory with external sources is 

exp{iZ(/^,r^, = j VqVqVG^ exp |i J d'^xCQCD{q,q,G^,lf„r^, s,p)^ (21) 

At low energies, i.e. below 1 GeV in QCD, the effective action can be obtained also 
from the low-energy effective theory 



exp{zZ(/^,r^, = y" VU exp |z J d'^xCeff{U,lf„r^,s,p)^ 



(22) 



With this effective action, the n-point Green functions can be easily derived by taking 
the functional derivative w.r.t. the external sources of Z{J) 



5" 



i'^6j{xi) . ..Sj{Xn) 



Z[J] 



J=0 



(23) 



Here j stands for any of the external sources l^,r^,s,p and J for the whole set of 
them. 



The vector current f ^ and axial- vector current included via 



ru = Vu + a'u- 



(24) 



In this paper we will calculate the two-point functions of vector, axial-vector, 
scalar and pseudo-scalar currents. The fermion currents in the complex case are 
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defined as 



V^{x) = q,^J,qj, (25) 

A;^{x) = q,Ttji,i,qj, (26) 

S'^ix) = -q,^q,, (27) 

P'^ix) = tq,^7,qj. (28) 



is an SU{n) generatoiH or in addition for tlie singlet scalar and pseudo-scalar 
current the unit matrix which we label by T°. These currents also exist the for 
real and pseudo-real case. In this case also currents with two fermions or two anti- 
fermions exist. These can be combined with those above. The generators can then 
become SU{2n) generators. All conserved generators are like the vector or scalar 
case while the broken generators are like the axial- vecor or pseudo-scalar case. All 
those cases are related to the ones with the currents of f l25p - (|25]) via transformations 
under the unbroken part of the symmetry group. 

The definitions of the two-point functions are 

IlvaM = 1 1 d'x e'"-^ {0\TiV^ix)V:iO)y\0) , 
UAa,Aq) = t J d'x e"^-"" {0\T{A;{x)At{0))^0) , 
IlsMaM = 1 1 d^x e^'^-- {0\T{V;ix)S''{0)y\0) , 

npMa^(g) = ^l^^xe^^-^ (o|r(A«(x)P'^(o))t|o), 

Usaiq) = tjd^xe'^-^ {0\TiS''ix)S''iO)y\0), 

Upaiq) = iyrfVe'^-" (0|r(P"(x)P"(0))t|0). (29) 

Using Lorentz invariance the two-point functions with vectors and axial- vectors can 
be decomposed in scalar functions 

^Vat^u = {qt^qu - q^gf,u)'n.'yl{q^) + qt^qu^vli<l'^) ■ (30) 

where ny^(g^) is the transverse part and n['?^(g^) is the longitudinal part or alter- 
natively the spin one and spin part. The same definition holds for the axial- vector 
two-point functions. The mixed functions can be decomposed as 

^SMafi = qtMsMa , 

HpMa/i = iqfJ^PMa ■ (31) 

Using the divergence of fermion currents and equal time commutation relations, we 
find that some two-point functions are related to each other by Ward identities. In 
the equal mass case considered here, they are 

nl/i = IIsAfa = , 
q^Yif^ = ImXlpua. 

q^V^fa = 4m2np, + 4m(gg). (32) 

^We have defined here the current with fcrmion-anti-fcrmion operators, hence the SU{n) for n 
fermions. For the real and pseudo-real case, the unbroken symmetry relates them also to difermion 
ot dianti-fermmion operators. 
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(4) (5) (6) 





(8) (9) (10) (11) 



"■^^^ ^^ ^^ ^■✓N^' 'N/N/^ ^^ ^^ ^^^ ^^ ^^ ^ ^^ ^-V/^ 



(12) (13) (14) (15) 

Figure 1: The diagrams for the vector two-point function. A filled circle is a vertex 
from £2, a filled square a vertex from £4, and an open square a vertex from Cq. The 
top line is order p^. The remaining ones are order p^. 



The vacuum expectation value is the single quark-anti-quark one. We will use the 
last relation to double check our results of axial-vector and pseudo-scalar two-point 
functions. 

The mixed two-point functions, HsMa and UpMa we do not discuss further since 
they are fully given by the Ward identities. 

3.2 The Vector Two-Point Function 



The vector two-point function is defined in fl29l) . The longitudinal part vanishes for 
all three cases because of the Ward identities. 

The Feynman diagrams for the vector two-point function are listed in Fig. [T] 
There is no diagram at lowest order. The diagrams at NLO are (1-3) in Fig. [H The 
NNLO diagrams are (4-15). The 3-fiavour QCD case is known to NNLO [221123]. 

We have rewritten the results in terms of the physical mass and decay constant. 
For these we use the notation Mm and Fm rather than the Mphys and -Fphys used in 
pT| [T2] . Their expression in terms of the lowest order quantities F and = 2Bom 
can be found in [TT]. We also use the quantities 

log — — and TTie = ■ (33) 



The loop integral B22 is defined in Appendix lA.ll 
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The results up to NNLO for three different cases are hsted below, where the first 
line in each case is the NLO contributions, the rest are NNLO contributions. 



n 



(1) 

vv 



Complex 



n 



AB22{Ml,Ml,q^) + 2LMl 



n 



(1) 

vv 



AM 



M 



M'^ 



L'n' - Sq'Kl,, + mi,{LL\^n - AK;, - AK^.n) 



Real 

q^ 



AB22{Ml,,Ml„q') + 2MlL 



AL[, - 8HI 



AM'^ 

^L{n + iy-W{n + l)L; 



B,2{Ml„Ml„q') 



+ -in + iy[B,2iMi,,Mt,,q')Y + 



q^ 



(34) 



Sq'Kl,, + 8Mij[LL[,{n + f ) - AK^, - 8K 



82 



n\ 



(35) 



n 



(1) 

vv 



Pseudo — Real 

-^{n - 1) [4522(M^, Ml, g2) + 2MI,L 
q2 L 

1 



AL[, - 8HI 



AM"^ 

^L{n-iy-16Ll{n-l] 



522(M^,M^,g^) 



+l(n - im^iMl, Ml, q-)f + ^L\n - If 



q^ 



-8q'K[,, + 8Ml[LL\,{n -I)- AK;, - 8Kl^n] 



(36) 



The complex result with n = 3 agrees with [22j when the masses there are set equal 



3.3 The Axial- Vector Two-Point Function 

The axial- vector two-point function is defined in (129|) . Similar to the vector two-point 
function, it also can be decomposed in a transverse and longitudinal part. 

n^A = {<i'<f - q'gn^^lW) + g'^g^nJl(g2) . (37) 

The diagrams contributing at LO are shown in (1-2) in Fig. |2l The LO results 
are the same for all three cases. The result is 

^'AA{q') = 1^9^'^ - q'q'-r^)j ■ m 

F and M are the LO decay constant and mass respectively. Note that in the massless 
limit this has only a transverse part as follows from the Ward identities. 

The diagrams at NLO are (3-10) in Fig. |2]and the NNLO diagrams are (11-48) 
in Fig. [31 
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(1) 



(2) 



^ — ■ — ^ 



(3) (4) (5) (6) 






(7) (8) (9) (10) 



Figure 2: The axial-vector two-point function at LO and NLO. The filled circle is a 
vertex from £2, The filled square is a vertex from £4, and the open square is a vertex 
from Cq. 



Many of the diagrams are one-particle-reducible and at first sight have double 
and triple poles at = M^. From general properties of field theory these should be 
resummable in into a single pole at the physical mass, = and a nonsingular 
part that only has cuts. The residue at the pole is the decay constant squared. We 
must thus find contributions that allow for the last term in fj38|l the lowest order 
F^, to be replaced by F^, Mf^. It turns out to be advantageous to also do this 
in the first term. Most of the corrections are already included in this way. 

At NLO the remaining part is only from the tree level diagram (3) in Fig. |5]and 

is 



n 



(1) 

AA 



AL[, - 8HI , n 



(0) 
AA 







(39) 



So we can express our result up to NNLO as 



Ka{<i' 




r(i) 



+ (g'^g^-gV'^)(4L^o-8i/[) 



(40) 



The n^]j(g^) and Ii^^j^{q^) are the remainders at NNLO and have no singularity at 
«2 - /If 2 

The transverse part can be obtained from the part containing g^^ as an overall 
factor. So the transverse part cannot come from the one-particle reducible diagrams 
and only gets contributions from diagrams (11-16) at NNLO. The sunset integrals 
and H21 appearing in the results are defined in Appendix IA.2I 
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(11) 



(16) 




(21) 




(26) 



(31) 



(36) 



(41) 



(46) 



(12) 



(17) 




(22) 




(27) 



(32) 



(37) 



(42) 



(13) 



(18) 




(23) 




(28) 



(33) 



(38) 



•N/N^^ 0^ 

(43) 





(14) 



(19) 




(24) 




(29) 



(34) 



(39) 



(44) 



(47) 



(48) 



(15) 




(20) 




(25) 




(30) 



(35) 







(40) 



— ■- 



(45) 



Figure 3: The axial-vector two-point function at NNLO. The filled circle is a vertex 
from C2, the filled square is a vertex from £4, and the open square is a vertex from 
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The longitudinal part gets at NNLO contributions from all diagrams shown in 
Fig 131 In order to rewrite the results into the single pole we need to expand the 
integrals around the mass. This introduces derivatives of the sunsetintegrals. These 
always show up in the combinations if*^ and defined in Appendix IA.2I 



Complex 



Y 



M 



^'H^iMl,, Ml,, Ml,, q') - H!,{Ml, M^, M^, 



2^2 



~6 



SnLLlo - 32 {Kl^^ + K{^^n + Kl, + K 



18 



n] 



-q\lQK{,, + m[ 



115> 



+ ^-7T + 



+71 



16 



Ml, 



q^ 

+mW 



7 1 9 9 I7n 

— n'^Ti'^ 

2^2 8 64 



6 
6 



-n 



IT 



Ml, 



V36 ' 72 
4 4 



n 

96' 



(41) 



2g2 



H'' {Ml, Ml,, Ml, e 



^Ml, 2 zjM( nj2 ji^2 Ji/f2 n'^\ j_ ^MIl T^-^ 
^21 l^^^M? M]^, M^, q ) + —^^113 

^2 , (7 15^2 



Real 

n 



/2 n2 2 

^^^^ U^y-;^ 



16 



6 



64 



2^2 



(42) 



(n+ 1) 



M 



*^ H'^iMl, Ml,, Ml, q') - H^,{Ml, Ml, Ml, q' 



M 



■ r2 

—n{n + 1) - 8nLL^o - 32(ir[o2 + 2nir[o3 + ^Ir + ^nK{^] 
6 



-g2(16i^[o9 + 8^ 



T N 

iisJ 



g2 



^K[i3 + -n(n + l)L' 



Ml, /I3n2 



+^?6 



Ml, 



1 43n 
+ — + 



+ 



1 2 2 

--n TT 



2n2 24 2n 6j ' 6 
2 17^2 7 5n 7 



nn 



64 ^ 8r22 192 8n 24 



+Mln{n + l) ^— + l^+^(n 



96 



i)g' 



(43) 



g2 2 



-|Ml,n(n + l)if2^ {Ml, Ml, Ml, q') + °-^K[,, 
z q 



H'' {Ml, Ml, Ml, q' 
8Mf 
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-TTiQ- 



M 

7 



M 



n 



L\ — 



_^^2 Ml, 

+-16^ 

Pseudo — Real 



2n2 
7 



7n 
24 



1 

2n 



6 



lOln 7 7 
192 8^ ^ 24 



(44) 



n 



[n 



M 



M 



H^{Ml, Ml, Ml, q') - Hl{Ml, M^, M^, q' 



-32Ml{K 



+ 2nir[o3 + K[, + 2nK 



8Mtr^^^ Ml3n 



-"■113 



n-l)L^ + 



is; 



Ml, 



q^ 
13^2 



6 



1)L^ - 8MlnLL\Q 



2^2 



43n 1 1 
^ ^ 2^ ~ 6 



M 



n 



n 
6 



-'-2 2 



+M^n(n - 1) 



TC 



17n^ 7 nn^ 

2 1 



5n 7 7 ' 
192 ^ 8^ ~ 24, 



36 72 



n 



+ ^96^" 



(45) 



M6 



M 



n 



2 2 



in-1] 



MtAn + l){ — -^ 



--Mln{n-l)H^i{Ml,Ml,Ml,q') + 
z q 



H'^' {Ml, Ml, Ml, q') 
8Mf 



■M T^r 
2 -"-113 



+vri6- 



M 



M 



L 



n 



1 7n 
2^^2A 



1. 1 

2^ ^ 6 



^2 Ml 



15^2 



7 lOln 



192 8^ ^24, 



(46) 



The axial two-point function is known in 3-flavour ChPT 



We have 



checked that our result agrees with the one in [22] in the limit of equal masses. 



3.4 The Scalar Two-Point Functions 

The scalar two-point function is defined in (12^ . which contains the unbroken gener- 
ator case {T"- = Q"-) and the singlet case (a = 0). 

The Feynman diagrams for both cases are the same as those for the vector two- 
point function shown in Figure [1] except that diagrams (2) and (5-7) are absent. 
Diagrams (1) and (3) are at NLO, and the diagrams (4) and (8-11) are at NNLO. 



3.4.1 case 

The scalar two-point functions are similar to the vector two-point functions, the LO 
results are zero for all the three cases since the vertex at LO is absent. We have 
rewritten again everything in terms of the physical mass and decay constant, Ml 
and Fm- The results for the three cases are given below. The first line is the NLO 
contribution and the remainder is the NNLO contribution. 
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Complex 



+MlL 



64 

n 



- 32n^ LI - 64L; + (n^ - a) ^Ll 



Real 



8g' 



32 



64 



n 



n 



-^Ll + Mi, I —L - 16L: - 32L^ + —L 

n 

(q^ 2M' 



2 ' 
M 



Uss = B',\8H; + 16Ll + -{n-l){n + 2)B{m\q') 



+MlL 



32 



n 



16 



32 - 32n - 64L^ + {n - 1) {n + 2) — L; 



n 



+B{m\q^) (n-l)(n + 2) 



32 ..rr 64^/ 

--L^ + 64L^ + -L^ 



V^5 + Ml,| 



1 1 

- + - 

n 



L - 32L" 



^ + M^. ( — - ^ 
2n 



+B{m^,q'^f{n-l) {n + 2) 

Pseudo — Real 

Uss = Bl\8Hl + lQLl + ^{n + l){n-2)B{m\q^)^ 

q^ (8ir[i3 + 32ir^7) + Ml (l92Kl^ + 128nK, 
'32 



-4. 



26 



16 

32 - 32n ) LI - 64L; + [n + I) {n - 2) —L 



n 



+B{m\q^) (n + l)(r2-2) 



32^^ 64^, 
— + 64L^ + —LI 



+B{m^,q^f [n + l) {n-2) 



- + 4 U - 32L^ 



1 1 
2n 



(47) 



(4J 



(49) 



The definition of the one-loop function B[m?, q^) can be found in Appendix lA.ll 
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3.4.2 Singlet ceise 

We have also calculated the singlet case. This is the quark-antiquark combination 
that shows up in the mass term. 

We write the expression up to NNLO as: 

Complex 

= BU8nH; + 32n^Ll + 16nLl + 2{n'^ -l)B{m^,q'^)\ 



+MltL (n^ - l) {32nLl + 32L^ - 64nL^ - 64L^) 



16q^ {nil + LI) 



+Ml, (^L + 64 {2LI + 2nLl - LI - nLl)j 



+B{m^,qY {n^ - l) [nq 
Real 



2 2M|,' 
n 



(50) 



= BU IQnHl + 128n^Ll + 32nL^ + 2(2n^ + n - l)B{m^, q^^ 



Bl 



16q^ (ni^[i3 + 4nKl^ + Sn^Kl^) + 384M^ (71/^2^5 + 2n^K;Q + 4n^/s:;7) 



+M^L (2n^ + n - 1) (64nL^ + 32L^ - 128nL^ - 64L^) 
+B{m'^, q^) {2n^ + n - l) 



16g2 (2nL^ + LI) 



+Ml, (^(^-2 + ^) L + 64(2L^ + 4nL^ - - 2nL^)) 



+fi(m^g^)' [2n^ + n- 1 
Pseudo — Real 



nq' + Ml,(l- 



(51) 



= S^|l6r^i/^ + 128r^^L^ + 32r^L^ + 2(2r^^-r^-l)S(m^?^)| 



■4. 



+MIjL (2n^ - n - 1) (64nL^ + 32L^ - 128nL^ - 64L^) 
+S(m^ g^) (271^ - n - 1) 



16^2 (^2nLl + LI) 



(^(^2 + L + 64(2L^ + AnL^ -LI- 2nL 
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— n 



nq^ + Ml 



4-1 

n 



(52) 



We also written the result in term of physical Mj^ and Fm- Notice that all loop 
diagrams are proportional to the number of Goldstone bosons in each case, i.e. n^ — 1, 
2n^ + n — 1, 2'n? — n — 1 for the complex, real and pseudo-real case respectively. 



3.5 The Pseudo-Scalar Two-Point Functions 

The pseudo-scalar two-point function is defined in fl29|) . Just as in the case of the 
axial-vector two-point function there are one-particle-reducible diagrams. The dia- 
grams are the same as those for the axial-vector two-point function with the axial- 
vector current replaced by a pseudo-scalar current. These are shown in Figure [2] and 
El There is also no vertex with two pseudo-scalar currents at LO so the equivalent of 
diagrams (1) and (7) in Figure [2] and (13-15) in Figure |3] vanish immediately. Just 
as in the scalar case, one should distinguish here between two cases: The adjoint 
case for the complex representation case which generalizes to the broken generators 
for the real and pseudo-real case, and the singlet operator with in ( 125]) the unit 
operator. 

In Section 13.31 we could simplify the final expressions very much by writing the 
final expression with the single pole at the meson mass in terms of the decay constant. 
The same happens here if we instead rewrite the result in terms of the meson pseudo- 
scalar decay constant Gm- So we first need to obtain that quantity to NNLO. 



3.5.1 The meson pseudo-scalar decay constant Gm 

The decay constant of the pseudoscalar density to the mesons, Gm is defineclfl simi- 
larly to Fm- 

m^.T'^qln') = -^5^'Gm (53) 

The calculation of Gm is very similar to Fm, the diagrams are exactly those shown 
in Figure 2 in [TT] with one of the legs replaced by the pseudo-scalar current. There 
is here also a contribution from wave-function renormalization. In [11] we reported 
all the quantities Mlj, Fm and (qq) as an expansion in the bare or lowest order 
quantities F and = 23^171. We therefore do the same here. We therefore use the 
quantity 

1 

instead of L as in the other sections of this paper. 

This quantity has been calculated to NLO in two-flavour ChPT in [IB] and was 
called Gt, there. We have checked that our NLO result agrees with theirs. 

At leading order, all the three case have same expression: 

^,,, = ^0 = 250^. (55) 



*The \/2 is included in the definition to have the same normahzation as [18]. 
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We express the full results up to NNLO in terms of the LO meson mass and 
decay constant F as 



Gm = 2SoF h + _ ac + 6g 
At NLO and NNLO, the coefficients aa and ha are 



(56) 



ho = 



ac = 
ha = 



Complex 

1 n 



n 



Lo + A{-nLl -LI + AnLl + 4L^) 



8K[, 



-64(L^ + nL^iLl + nil) + 24(L^ + nL^)' 



'25 



39 



^18 



+7ri6 



- (32 - 22n2) (^Ll + ^L^) + (4 - 8n^){Ll + AL] 

(8 + 2n 

(2 - n') (-LI + 8L^6 - -LI - - -L 
\n n n 



80 



48n L: + 



(- 



lOn L 



12 



4n L: 



+Ll 



_13 
V 256 ' 24 ^ ' 
^3n2_3 _9_^ 
^2^ 



13 



55n 
~96" 



- 1 + 



2n2 



Real 



16 
1 

2 + 2 



(57) 



1 

2n 



- 4LI + 32nLl + 16L^) 



64(L^ + 2nLl){Ll + 2nL^) + 24(L'5 + 2nL^)2 

k;, + 48K;, + 32K;, 



-32Kl^n^ 



(-16 + 16n + 2271^) (2LI + ^Ll^ + (4 - 877 - 16n'^)Ll 
+16(1 - 377 - 4772)L^ - (^40 - ^ + 4877^ - (^6 - ^ + IO77) L 
- (8 + 277 + 477^) LI - (6-^ + 477) LI 

'l-n-n^)( -L: + 16L^ - -LI - 8L: - -LI 
\n n n 



+(77 + 2772)L^ + 2L\ 

/II3772 44377 13 
256 



+ 277 Un 



77 



+ 



13 



768 9 6 3277 



13 \ 
32^j 
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be 



55rr 



96 

^9 / 3n^ 13n 

iTiT + IF 

Pseudo — Real 



Q7n 3 5 
"96" ~ 8 ~ 8^^ 



1 



3 9 

+ 



n 



1 



1 



8 2n 8n2 



:-8nL: - 4Lr + 32nL" + IQLl 



r\2 



(58) 



-64(L^ + 2nL^)(L^ + 2nL^) + 2A{Ll + 2nLl) 



32K';, - 8Klg - 8K;^ + A8K 



23 



'■25 



+32^3^9 



17 

16 - 16n + 22n^) (2LI + -Ll) + (4 + 8n - 16^^)^^ 



n 



+ 16(l + 3n-4n^)L^+ 40 + 



40 



n 



48n + 6 + 



6 



n 



10n]L 



- (8 -2n + 4n2) + ( 6 + - - 4n ) 



6 

n 



+vri6 



[l + n 



n 



8l: 



-a 



n 



n 



n 



+ {-n + 2n^)L^ + 2L^ + ( -1 - - + 2n ) L; 



113^2 443n 



256 
/ 55n^ 



768 
67n 



13 

96 ^ 32n 



13 13 \ 
^ 32^J 



^ 96 96 
3n2 13r2 1 



3 

- 77 + 



5 7 



5n 



16 



16 8 



(59) 



3.5.2 X" case 

The pseudo-scale two point functions are similar to the axial-vector ones in the 
diagrams as described above. The LO result is the same for all the three cases: 



TTO 
iipp 



ri2 



2 g2 _ M2 ■ 



(60) 



The superscript "a" indicates the case with T" in f l28|) an SU {n) generator. For the 
real and pseudo-real case this is related by the conserved part of the symmetry group 
also to a number of diquark currents. 

Subtracting the pole contribution in terms of the physical mass and decay con- 
stants, Mf^, Fm and Gm, absorbs the major part of the higher order corrections. 
The final results are thus much simpler when written in this way. The remaining 
part at NLO is 

16L^) . (61) 



TTO 
iipp 



Bq(8HI 
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Thus we can define the full NNLO results as 



Bl 



(62) 



where the Ilpp is the remainder at NNLO. Its expression for the three different cases 
is: 



LLpp 



Complex 

-\n^q'H^,{Ml,,Ml,Ml,q-) 



+ 



n 



Q 



n 



Real 



3 12 
5 _ 85n2 
3 



H^{Ml,,Ml,Ml,q^) 
+ 2 ) + LMIj (ULI + 32nL^ - 



n 



96 



7 15n^ 



64 



2^ 



(63) 



rra 
LLpp 



q%{n + l)H^I{Mi,,Mi„Mi,,q') 



rra 
LLpp 



+ 



1 n 1 1 

+ - + 



3 n 3 
+8q'Kl,, + 64Ml,{K[, 
2 3 



q^- M'^q^n{n+1) 



K. 



39 



2nirio) 



+7^161/ 



n 



V 2 2n^ n+|^ + i)+32MiL|2LJ + L;„ 

' . / 2 7n 2 2\ „ / 



+7^16 



Ml 



Pseudo 



12 

(_S5rf_ _ 
\ 96~ ~ 

15n^ _ _]_ 
8n2 



12 

5 _ 125n 
4^ ~ ~96~ 



n 



1 7n 1 r 
2^~ 24 ^2^^ 6, 



^ 4n ^ 12 . 



64 
Real 

rM/ 



lOln ^ 
192 ^ 8n ^ 24, 



(64) 



^qMl-n)H^,{Ml,Ml,Ml,q') 



+ 



1 n 1 1\a 1,^9 9/-, X 



{Ml, Ml, Ml, q') 



+Sq'Kl,, + 64M^(ir[7 + 2nir[8 - i^g^g - 2nKl,) 
' 3 3 1 



+MIL 



2n2 + 2n + 2 



+ MlL 



+ 32 n - 



n 



l]Ll 



, . / r?/ 2 7n 2 2 \ o 1 7n 1 1 



12 



12 



8 2n2 24 2n 6, 
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.2 / 85^2 5 125n 5 5 
^ I I 

^ \ 96 4n2 96 An 12 
15n2 7 lOln 7 7 \ 
"gT ~ 192 ^ 8^ ^ 24 j ^ 



(65) 

The loop integrals if and i/|f are defined in Appendix IA.2I 



3.5.3 Singlet case 

In the singlet case with a = 0, there is no contribution with poles. Only the one- 
particle-irreducible diagrams contribute. As a consequence, there is no order p"^ 
contribution and at order p'^ there is only a tree level contribution from the equivalent 
of diagram (3) in Figure [2J At order p^ or NNLO only the one-particle-irreducible 
diagrams contribute and since there is no order vertex with two pseudo-scalar 
currents only diagram (11-12) and (16) in Figure |3] contribute. 

Since there is no single pole contribution, there is also no need here to expand in 
the integrals around the meson mass. The integral is defined in Appendix IA.2I 

The singlet pseudo-scalar two-point function we write as 

U'pp = B^Tfpp + ^fl° p . (66) 
The results for the three cases are 



Complex : 

iH 
2 



+L^Ml^- (n^ - l) (n^ - a) + 6ALMlj{n^ - l)(nL^ + L^) 



-FM>?,l(n^-l)(n^- 4) ^ + 1 , (67) 



Real : 

2 



Hpp = 16ni72^ - 128^2^^ - 32nL^ , 



11% = -—(2n' + n-l)[n^ + n-2)H^{Ml,Ml,Mlj,q^) 

+g2 (iQKl-^^n - 128Kiy) - 128M^^ (hK;^ + 2n^Kl^ + 2n^Kl^ + An^Kl^ 
+L^MIj- (2n^ + n - l) (n^ + n - 2) + 6AMl^L{2n^ + n- l){2n^ + L^) 



n 

+M|^7r2gl {2n^ + ^ - l) (n^ + n - 2) + 1 j , (68) 

Pseudo — Real : 

TTpp = IQnHl - USn^Uj - 32nL^ , 

flip = -^{2n'-n-l){n'-n-2)H^{Ml„Ml„Ml„q') 
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Figure 4: The one-loop oblique correction to LEP process e"*" + e — )■ g + g. 

+L^MIj^ {2n^ - n - l) (n^ - n - 2) + 64:MljL{2n^ - n - l)(2nL; + L^) 
+Ml7rl^ {2n^ _ ^ _ l) (^2 _ ^ _ 2) + 1^ . (69) 

Notice that just as for the scalar singlet two-point function, all loop contributions 
are proportional to the number of Goldstone bosons. 

3.6 Large n 

As one can see from all the explicit formulas, many of the expressions become equal 
for the different cases in the large n limit . 

4 The Oblique Corrections and S-parameter 

The physical process at the CERN LEP collider is + e~ ^ q + q. There are 
three types of one loop correction to this process: vacuum polarization corrections, 
vertex corrections, and box corrections. The vacuum polarization contribution is 
independent of the external fermions and it dominates the contributions from physics 
beyond SM. For the light fermions, the other two corrections are suppressed by 
factor of rnj/rn^. That's why the vacuum polarization corrections are called "oblique 
corrections,", and the vertex and box corrections are called "nonoblique corrections." 

The oblique polarization only affect the gauge bosons propagators and their mix- 
ing. The vacuum polarization amplitude can be defined as 

g^'UxY + (gV terms) = i j <fx e'^'" (0|T(Jj^(a;) J{;(0))|0) . (70) 

The influence of new physics to the oblique corrections can be summarized to three 
parameters: 5, T and f/. One can find their definition in Ref. [15]. These parameters 
are chosen to be zero at a reference point in the SM. In the past 20 years, they have 
been studied intensively in many models beyond the Standard Model physics. 

For a beyond the Standard Model with strong dynamics at the TeV scale, there 
will in general be many resonances and other nonperturbative effects. At low mo- 
menta one can use the EFT as described above for these cases. In this paper, we will 
estimate the S parameter contribution from pseudo-Goldstone Boson sector within 
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the EFT. The parameter T and U vanish because of the exact flavor symmetry, i.e. 
we work in the equal mass case. 

The S parameter can be written a^ [T5] 



S = -2n 



nW(0) - n:,^(0)] = 2vr A (g^nW _ q^U^j\)^^^^ . (71) 



nyy(O) and n^^(O) are the derivatives of the vector and axial- vector two-point func- 
tions at = 0. One should keep in mind that S is defined to be zero at a particular 
place in the standard model, as discussed at the end of section V in [15]. Our formulas 
are the equivalent of (5.12) in that reference. 
The result can be written as 

S = S+^S, (72) 



with 



Complex : 



3 

S = U{Kl,,-Kl, + Kl, + nKl,,-nfKl^ + nKl,) + -L^ 



+16n {LI + 2L\,) L - ^i,^L + vr?y (— - -^J (73) 
Real : 

S = QA {Kl,, - Kl, + Kl, + 2nK[o3 - 2nfKl, + 2nKl,) + Ilil^L^ 
+16 [(n + l)Ll + {2n + l)^;^] L - ^^ /M^+l) ^, 

+->(«+l)(^-^^). (74) 
Pseudo — real : 

T7 . ^ r r 2(n — 1)tT , ^ . 

S = -levrLio - ^ g ' {L + TTie) , 

^ = 64 (ir[o2 - K;, + Kl, + 2nK{,., - 2nfK;, + 2nK{,) + ^ ^ ' 

+16 [{n - l)Ll + {2n - l)L^o] L - vne ^ 

-l)(^-|^). (75) 

The quantity if) is 

7/, = 6^3012 f ^) = 7.0317217160684 . (76) 



^Our two point functions are normalized differently from those in |15] . 
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Figure 5: The S'-parameter for the values of L'g and L'[q given in the text for the 
complex case, (a) n = 2 (b) n = 4. 



n=2 real 



n=4 real 
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Figure 6: The S'-parameter for the values of Lg and L^q given in the text for the real 
case, (a) n = 2 (b) n = 4. 



The real purpose of fl75]) - fl75]) is to be able to study the S'-parameter in more gen- 
eral theories than just scaling up from QCD. However to provide some feeling about 
numerical results we choose parameters as if they are scaled up from QCD/ChPT. 
We change = 0.0922 MeV to Fm = 243 GeV and the subtraction scale from 
0.77 GeV to 2 TeV. We set the = and keep LI = 0.00593 and L^q = -0.00406 
at their values from ChPT [2S1 121]. 

In Figures O |6] and [7] we have shown the results for our three cases complex, real 
and pseudo-real for n = 2 and n = 4. Shown are the full and contributions 
as well as the part proportional to L^q only. The latter is what is the usual 
contribution to S corrected for the pieces that go into the reference point at p^. We 
cannot do the same for the full result since that depends on how one treats the extra 
pseudo-Goldstone bosons that occur in the other models. 
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Figure 7: The S'-parameter for the values of Lg and L'Jq given in the text for the 
pseudo-real case, (a) n = 2 (b) n = 4. 



5 Conclusion 

In this paper, we have calculated the two-point correlation functions of vector, axial- 
vector, scalar and pseudo-scalar currents for QCD-like theories. 

In the beginning of the paper, we gave a very brief overview of the QCD-like 
theories and their EFT treatment as developed earlier. 

We then gave the analytic results of those two-point functions up to NNLO. The 
results are significantly shortened by using the physical meson mass M^j and decay 
constants Fm and Gm when rewriting the pole contributions. 

The main use of these formulas is expected to be in extrapolations to zero fermion 
mass of technicolour related lattice calculations. We have therefore also included 
precisely the combination needed for the S'-parameter. 
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A Loop integrals 



We use dimensional regularization and MS scheme to evaluate the loop integrals, 
ci = 4 - 2e. 
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A.l One- loop integrals 

The loop integral with one propagator is 
2^ _ 1 f d^q 1 



m 



167r2 ^ 



Xq — In(m^) + e 



Here 



C = ln(47r) + 1-7 



Ao = - + C 
e 



The extra +1 in C is the ChPT version of MS. 
The loop integrals with two propagators are 



TDt 2 2 2\ 

Bf'imlmlp) 



iJ {2nY {q'^ — 'mi){{q — py — ml) 



2\ ' 



d'^q 



{27ryiq^-ml){{q-py-ml) 
p'^Biimlmlp'') , 

d'^q q^'q" 
{27rr{q2-ml){iq-py-ml) 



1 r d'^c 
1 J (27r 



= p^p^B^iimj, ml,p^) + g^''B22{ml, ml,p^) . 
The two last integrals can be reduced to simpler integrals A and B via 

Si(m^m^p2) = ^B{ml,ml,p'^) , 



TD / 2 2 2\ 

B22[m ,m ,p ^ 



2(rf- 1) 



m^) + 2m^ - -p^ B{m% m,p^) 



B. 
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(77) 



(78) 



(m^m^p^) = ^ [A(m^)+m^S(m^m^p^) -dS22(m^m^p^)] . (79) 



167r2 
1 



Ao + S(m^p^) + C»(e), 



We quote here only the equal mass case results relevant for this paper. The explicit 
expression for B is 

B[m ,m ,p ) — 
B{m^,p^) 

J{m'^,p'^) 
The function J{m'^,p^) is 



167r2 

1 Y 



-1-m log — + J(m2,p ) , 



167r2 



dxln 



' — x{l — x)p^ 



(80) 



J{m'^,p^) 



2 + aln(f^), p^<0, 

2 - 2^1 - 1 ■ arccot - l) , < < irn^^ 

2 + (T In (^i^) + ivra, > Am"^, 
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/ 4 

a(x) = ^l--, x = — ^[0,A]. (81) 

Taking derivatives w.r.t. at = is most easily done in the form with the 
Feynman parameter integration exphcit. 

A. 2 Sunset integrals 

The sunset integrals are done with the methods of [22| 128] . They are defined as 

1 f d'^q d'^r X 
^^^^^ ^1? J {2tiY {2t:Y (g2 - ml) {r^ - ml) [{q + r- pf - ml] ' ^^^^ 

The various sunset integrals with Lorenz indices are 
H{mlmlmiy) = ((1)), 

H^'imlmlmlp') = ((g^)) = p^/7i(m?, m^, m^; p^) , (83) 
H^''{mlmlmlp') = {{q^q^)) 

= p^'p''H2i{ml,ml,ml;p'^) + g>"'H22{ml,mlml;p'^) . 

and 

((r^^)) = p>'Hi{ml,ml,ml;p^) , 
{{r'^r")) = p^p''H2i{mlml,ml;p^)+g^^''H22{ml,ml,ml;p'), 

((gMr'^)) = p^^p'^H23{mlml,ml;p')+g^^H2,{ml,ml,ml;p'), (84) 

The function H is fully symmetric in m\,m\ and m\, while ifi, H21 and H22 are 
symmetric under the interchange of m\ and m\. The relation between the above 3 
functions 

2tj f 2 2 2 2\ I jTj I 2 2 2 2\ 

p H2i{mi,m2,m^;p ) + dH22{mi,m2,m^]p ) = 

m\H{m\, ml, ml; p^) + A{ml)A{ml) , (85) 

allows to express H22 in terms of H2i- 

Similar to the integral B and there is also a relation between H and Hi which 
in the equal mass case becomes 

Hi{m\ m\ m^- p^) = \H{m\m\ m"- p") . (86) 

The other functions, H23 and H23, can be written in term of H, Hi and if 21 by using 
relations derived from redefining the momenta and masses in its definition [22j. 

The full sunset integral expressions and the definition for finite part Hf = 
{H^,H[,H2i} can be found in the appendix of [22]. In our case we take mi = 

7712 = fTT-S = 
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In order to eliminate the extra poles in the expressions, sometimes we need to 
expand the Hf (m^, m^, m^; q^) around the pseudoscalar mass m^, and we define 



ttM/ 2 2 2 2\ 

ii j [m ,m ,m ;q ) 



ttF/ 2 2 2 2\ ttF/ 2 2 2 2\ 

ii j [m ,m ,m ; q ) — [m ,m ,m ■,m ) 

(87) 



/ 2 2\ TjFlf 2 2 2 2\ 

—(5 —m [m ,m ,m ■,m ) 



where 



H[ (to , m , m ; m ) = 



dq 
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